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SUMMARY 


Unsteady rotor wake interactions with the empennage, tail boom, and other 
aerodynamic surfaces of a helicopter have a significant influence on its aerodynamic 
performance, the ride quality, and vibration. A numerical method for computing the 
aerodynamic interaction between an interacting vortex wake and the viscous flow 
about arbitrary two-dimensional bodies has been developed to address this helicopter 
problem. The method solves for the flow field velocities on a body-fitted computa- 
tional mesh using finite-difference techniques. The interacting vortex wake is 
represented by an array of discrete vortices which, in turn, are represented by a 
finite-core model. The evolution of the interacting vortex wake is calculated by. 
Lagrangian techniques. The viscous flow field of the two-dimensional body is calcu- 
lated on an Euler ian grid. The flow around circular and elliptic cylinders in the 
absence of an interacting vortex wake has been calculated. These results compare 
very well with other numerical results and with results obtained from experiment and 
thereby demonstrate the accuracy of the viscous solution. The interaction of a 
rotor wake with the flow about a 4:1 elliptic cylinder at 45° incidence was calcu- 
lated for a Reynolds number of 3000. The results demonstrate the significant varia- 
tions in the lift and drag on the elliptic cylinder in the presence of the interact- 
ing rotor wake. 


SYMBOLS 


C d drag coefficient 

C fx force coefficient in x-direction 

Cfy force coefficient in y-direction 

lift coefficient 
J transformation Jacobian 

«, body length 

P transformation coordinate control parameter 

p pressure 

Q transformation coordinate control parameter 
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r radial variable 

r c vortex core radius 

r Q vortex core outer radius 

Re Reynolds number, U^fc/v 

t time 

u velocity component in the x-direction 

v velocity component in the y-direction 

x Cartesian coordinate 

y Cartesian coordinate 

a coordinate transformation coefficient 

8 coordinate transformation coefficient 

Y coordinate transformation coefficient 

r circulation 

n coordinate in transformed space 

9 incidence angle 

v kinematic viscosity 

5 coordinate in transformed space 

4> azimuth angle 

u> vorticity 

Subscripts 

( ) interacting vortex wake quantity, derivative with respect to 

( ) reference value 

( ) 1 value at body surface 
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Superscripts 


( )' dimensional quantity 


INTRODUCTION 


The numerical prediction for the strong interaction between vortical wakes and 
the viscous flow field about bodies is of considerable importance in the design 
analysis of rotorcraft. The flow field surrounding a helicopter configuration is 
highly complex because of the lack of symmetry and the unsteady aspects of the 
flow. This unsteady, asymmetric flow field is complicated further in that the shed 
wake from the rotor blades interacts with other components of the aircraft (e.g., 
the rotor blades, fuselage, tail boom, and tail rotor). In general, these are 
strong interactions in which the rotor wake flows onto or passes very close to the 
other components of the aircraft. Vortex interactions occur at many flight attitudes 
(i.e., hover, descent, and low speed flight) and therefore have profound effects on 
the overall aerodynamic efficiency of the rotorcraft. The effects of the vortex 
interactions are realized as increased vibratory loading on the fuselage, decreased 
payload capabilities, and increased noise. 

Because of the complexity of the problem and the limitations of the computer in 
modeling flow field, the research effort addressing vortex interaction for rotor- 
craft has developed in two directions: blade-vortex interaction and rotor-fuselage 

interaction. Most of the vortex interaction research for rotorcraft has been 
directed toward blade-vortex interaction in forward flight. Blade-vortex research 
has concentrated on the flow field near the rotor tip on the advancing side of the 
rotor. The flow field in this region is transonic and the rotor blade is at near- 
zero incidence. Navier-Stokes calculations of blade-vortex interaction are being 
conducted with increasing efficiency (refs. 1 and 2). These calculations model 
limited separation because of the incidence of the blade and flow velocity. In 
addition, the blade-vortex interaction is generally not a direct interaction between 
the rotor blade and the rotor wake, but a very close passage. The exception is 
descent in which the rotor wake is convected through the rotor and directly inter- 
acts with the rotor blades. The flow beneath the rotor near the fuselage is quite 
different. Collisions between the rotor wake and the fuselage are present every- 
where, the flow is subsonic, and regions of large separation exist. 

In contrast to blade-vortex interactions, which occur near the rotor tip, 
vortex interactions beneath the rotor occur on the entire rotorcraft fuselage. The 
rotor-fuselage interaction has received far less attention than the blade-vortex 
interaction. An overview of the rotor-fuselage interactional aerodynamics can be 
found in reference 3. These interactions, which exist for most flight configura- 
tions, are caused by the impingement of the convected rotor wake onto the fuselage 
components. In contrast to the blade-vortex interaction, the flow field of the 
rotor-fuselage interaction is subsonic, highly separated, and vortex-dominated. The 
fuselage is immersed in the rotor downwash, which is an unsteady flow with very 
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large velocity gradients just inboard of the rotor tip. The fuselage is not stream- 
lined in the direction of the downwash; as a consequence, separated wakes form at 
the surface. These separated wakes are very unsteady because of the interaction of 
the impinging rotor wake. The result of the vortex interaction is the coupling of 
the vorticity produced at the surface of the fuselage with the interacting vortex 
wake. The flow becomes increasingly unsteady owing to the mutual convection between 
the interacting vortex wake and the surface vorticity. The resulting flow field is 
very complex, and therefore difficult to model and predict. 

The analysis of rotor-fuselage interactions as well as blade-vortex interac- 
tions is complicated further by the computational restriction on mesh size. The 
effects of the mesh size appear in the available computer memory, the computational 
effort, and the ability of the numerical method to accurately describe the interact- 
ing vortex wake over long periods of time. Currently, the available computer memory 
is inadequate for the three-dimensional mesh necessary to model the flow field 
surrounding a representative rotorcraft. Therefore, the computational mesh enclos- 
ing a body is generally concentrated near the body to resolve the flow at that 
location and is considerably less dense away from the body. Owing to numerical 
diffusion, the coarse mesh away from the body is inadequate for accurate predictions 
of the vortex dynamics surrounding rotorcraft. To calculate the rotor-fuselage 
interaction, a method is needed to accurately predict both the vortex dynamics away 
from the body on a coarse mesh and the separated viscous flow near the body surface. 

In this paper, a method is presented to calculate the unsteady, two- 
dimensional, incompressible interaction of a vortex wake with the separated flow 
about bodies. The method solves the velocity/vorticity formulation of the Navier- 
Stokes equations. The computational mesh is concentrated near the body surface to 
resolve the boundary layer and is increasingly coarse further from the body. The 
interacting vortex wake is modeled as an array of finite-core vortices (refs. 4 
and 5). The core radius is variable and is independent of the mesh spacing. The 
finite-core model eliminates the numerical diffusion associated with the coarse mesh 
spacing and provides for the accurate prediction of the vortex wake dynamics away 
from the body, where the flow field is inviscid and dominated by the rotor wake. 
Closer to the body, the flow is viscous and the convection of the rotor wake can no 
longer be considered inviscid. In this region the rotor wake interacts with the’ 
viscous separated wake and can no longer be represented by a finite-core model. The 
vortex core is acted upon by shearing forces which will contort and possibly sepa- 
rate the vortex core into several separate vortices. At this point the finite-core 
vortex is distributed to the computational mesh and allowed to convect as part of 
the viscous solution. 

The solution is started impulsively at t = 0. The interacting vortex wake is 
inserted into the flow field at a selected time and position. This method has been 
applied to the solution of flow about circular and elliptic cylinders without the 
presence of an interacting vortex wake and for the interaction of a vortex wake with 
the flow about an elliptic cylinder at 45° incidence. 
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Problem Formulation 


The vortex interaction problem is modeled by the velocity /vorticity form of the 
unsteady incompressible Navier-Stokes equations. The governing equations are 

for the velocity, 


V u = -<o. 


V 2 v = 


<D 


(1a,b) 


and for the vorticity, 


id. + (u<d) + ( voj) = 1/Re(V 2 w) (2) 

u x y 

with 

2 nt .2 .. 2 . .2 2 
V = 3( ) / 3x + 3( ) /3y , 

( ) t = 3( )/3t, ( ) x = 3( )/3X, ( ) = 3( )/3y, 

( ) xx = 3( ) 2 /3 x 2 ( ) yy = 3 2 ( )/3x 2 


where (x,y) are the Cartesian coordinates, Re is the Reynolds number, and t is 
the time. The variables (u,v) are the Cartesian components of the velocity and the 
vorticity id is defined by 


The nondimens ional variables are written 

x = x'/a, y = y ' /a , u = u'/U^, v = v'/U^ 
o> = w'/(u /a), t = t'/(a/u ), Re = u a/v 

CO 7 ao 7 GO 


where 

a = the body length 

Re = Reynolds number based on body length 
U,,, = free stream velocity 
v = kinematic viscosity 

A body-fitted computational mesh is used in the finite-difference solution to 
the governing equations. The methodology developed in reference 6 is used for the 
mesh generation. The method solves for the curvilinear coordinates as solutions of 
an elliptic differential system. The nomenclature is that of reference 6. The 
transformed coordinates are solutions of the equations 
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C + C = PU.n) 
xx yy 


(4a) 


n xx + n yy = (4b) 

where (C,n) are the transformed coordinates, P and Q are functions used for coordi- 
nate control. 

In the transformed plane, equations (4a) and (4b) become 

“ x t5 ' 2Sx 5n *■ r \n * j2(Px 5 * V = 0 <5a) 

“ y C5 ' 2By Sn * Ty nn * j2(Py 5 * V = 0 (5b) 

The boundary conditions for equations (5a) and (5b) are the specified Cartesian 
coordinates at the body and the far-field boundary. The computational mesh and 
transformed coordinates are shown in figure 1 . 

The governing equations can now be written for the transformed plane. For the 
velocity components 



- 26u 

+ yU 

+ J 2 (Pu r 

+ Qu ) 

= — J ( X W - X W ) 

(6a) 

ZZ 

5n 

nn 

5 

n 

5 n n 5 


ZZ 

- 20V 

5n 

+ Y V 

nn 

+ J 2 (Pv^ 

. Q»„) 

= J(y w - y^w ) 
5 J 5 n 

(6b) 


and for the vorticity 

w + [y (uw) - y_(uw) l/J + [x (vu>) - x (vw) ]/J 

z n 5 5 n £ n n S 

= (aw - 20w + yui )/(ReJ 2 ) + (Pw_ + Qw )/Re (7) 

5n nn 5 n 

where 


( ) h = 3( )/3t, ( ) = 3( )/3 5, ( ) = 3( )/3n, 

( )^ = 3 2 ( )/35 2 ( ) nn = 3 2 ( )/3n 2 

( ). = 3 2 ( )/353n 

5n 

The boundary conditions for equations (6) and (7) are 
u = 0, v = 0 at the body surface, 
u = cos(e), v = sin(e) at the far-field boundary 


( 8 ) 
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w = 0 at the far-field boundary 


( 9 ) 


where 0 is the incidence of the body to the free stream flow. The solution is 
started impulsively at t = 0. Instantaneously after the impulsive start, t = 0+, 
the flow is defined by a combination of potential flow and a sheet of vorticity at 
the body surface. 

A third-order-polynomial distribution of vorticity about the vortex position is 
used to represent the vortex core (ref. 4). This distribution of vorticity is 

u>„ u (r) = r/(2*r 2 )[20/3(r 3 /r 3 ) - 10(r 2 /r 2 ) + 10/3] (10) 

vw o o o 

where u vw represents the vorticity of the interacting vortex wake and r Q is the 
outer radius of the vortex distribution. The core radius, the radial position of 
maximum tangential velocity, is equal to 0.6376 r Q . The influence of the inter- 
acting vortical wake is included in equations (6a) and (6b) through the x and y 
components of -(7 x u yw ) • The interacting vortex wake is modeled as an arr$y of 
finite-core vortices prior to the interaction. However, when the interacting vortex 
wake is less than a "prescribed" distance from the body surface, the finite-core 
model is no longer adequate to model the wake. At this point, the interacting 
vortex wake is distributed to the computational mesh and included as part of the 
viscous solution. 

The surface pressure is obtained by the line integral of the Navier-Stokes 
equation along the body. At the surface 


or 


p = -(Bu) - ao> )/(ReJ) (11) 

n n 5 


Pi = P, 



aaj^)/J]dn 


( 12 ) 


where p^ is the pressure on the surface and p Q is the reference pressure on the 
surface. 


The reference pressure is obtained by the line integral of the Navier-Stokes 
equation from the far-field boundary to the body surface along a constant n line. 

P 0 = J* max [(X U + y^v) t - (x^vw - y^uoj) + (yu^ - Bw^ ) /.( Re J ) ]dC + 1/2 (13) 

The force components acting on the body are calculated from the integration of 
the pressure and shearing forces at the body surface. The force coefficients are 



(14a) 


n 2 

(w.x )dn 
i n 

1 

(u^y^dn (14b) 

where u>^ is the vorticity at the body. 

The lift and drag coefficients are 

C & = -C fx [sin(e)] + C fy [cos(0)] 

C d = C fx [cos(e)] + C fy [sin(0)] 

where 0 is the incidence of the body with respect to the free-stream flow. 

Numerical Method 

The governing equations are solved by standard finite-difference methods on the 
body-fitted computational mesh. The mesh is concentrated near the body to resolve 
the flow at that location. 

The vorticity transport equation, equation (7), is solved using an alternating 
direction implicit method (ADI) (ref. 7). The ADI method solves the time-dependent 
vorticity transport equation by separating the equation into two implicit equations 
for a single time-step. 

The finite difference expression for the vorticity transport equation is 
(o> n+1/2 - a, n )/(At/2) + [y n (uu.)£ +1/2 - y 5 (uu))J]/J + [^(vw)" - x n (va>)" +1/2 ]/J 

= (aa >"* 1/2 - 20a>" + Y 0 > n )/(ReJ 2 ) + (Pco" +1/2 + Qu) n )/Re (16a) 

55 5n nn 5 n 

for the ^-derivatives, and 

/ n+1 n+ 1/2 \ / /»i_ / o\ r.r / , \ ^ 1/2 

(to - o) )/(At/2) + [y^uto)^ 

- y (uto) n+1 ]/J + [x (vto) n+1 - x (vo))" +1/2 ]/J 
5 n 5 n n 5 

= (o< o__ - 20to_ + yw )/(ReJ ) + (Pu)_ + Qoj )/Re (16b) 

55 5n nn 5 n 

for the n-derivatives. 


(15a) 

(15b) 


'fx 


(p-jy^ ) dri - (1/Re) 


C fy = ( P 1 x^ )dn - (1/Re) 
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The vorticity at the body surface in equation ( 1 6b) is calculated from equa- 
tion (9) whereas the vorticity at the body surface at the half time-step is calcu- 
lated from equations (16a) and ( 1 6b ) as 


n+1/2 

o> 


(1/2)(u> n + (o n+1 ) + At/(4ReJ 2 )[(Yw n 

nn 


^ ) ] + ( AtQ)/(4Re) [</j 


n+ 1 

"n 1 


(17) 


The convective velocity at the half time-step in equations ( 1 6a ) and ( 1 6b ) is calcu- 
lated from the average of the n and n + 1 values of the velocity. 

The velocity equations (6a) and (6b) are cast in a time-dependent form and 
solved by the ADI method. 

The influence of the interacting vortex wake is included in the velocity equa- 
tions through the negative curl of the vorticity distribution in the vortex core, 

-(V x u ). The vorticity distribution in the vortex core is defined in equa- 
tion (iW. Outside the outer radius, the vorticity is zero. 

The solution to the time-dependent problem is calculated by solving the vortic- 
ity transport equation, updating the interacting-vortex-wake position, and then 
solving the velocity equations. 


Computed Results 

Results for the viscous flow about two-dimensional bodies have been calculated 
with and without an interacting vortex wake. The flow around circular and elliptic 
cylinders in the absence of an interacting vortex wake has been calculated in order 
to compare results of the current method with those of other numerical studies and 
experimental investigations. The flow around the circular cylinder was calculated 
for a Reynolds number of 3000. The flow around the elliptic cylinder was calculated 
for Reynolds numbers 200 and 3000 at 45° incidence. 

The interaction of a vortex wake with the flow around an elliptic cylinder at 
45° incidence and a Reynolds number of 3000 was calculated. The interacting vortex 
wake models the periodic passage of a vortex generated by a rotor blade. 


No Vortex Interaction Case 

Circular Cylinder - The flow about a circular cylinder was calculated at a 
Reynolds number of 3000. This condition was selected because there exists data at 
this Reynolds number for comparison. In addition, the emergence of several second- 
ary vortices behind the cylinder has been calculated numerically and observed exper- 
imentally. The computational mesh is 81x241 with the far field boundary at 
20 radii. The variation of the drag coefficient with time is shown in figure 2. 

The time, t = 1.0, corresponds to the time required to travel a distance equal to 
the cylinder radius at a velocity equal to that of the free-stream velocity. The 
effect of the secondary vortices on the drag coefficient can be seen as the small 
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variations in the curve. The distribution of the boundary vorticity at the cylinder 
surface is shown in figure 3. These results are in excellent agreement with the 
results of the numerical investigation of reference 8. The distribution of the 
u-component of velocity along the x-axis behind the cylinder is shown in fig- 
ure 4. Experimental results from reference 9 and the numerical results of refer- 
ence 8 are included in the figure for comparison. The comparison in the data is 
shown to be quite good. The streamline contours for the developing flow around the 
circular cylinder are shown in figure 5. The flow patterns demonstrate the emer- 
gence of the secondary vortices and are similar to those reported in references 8 
and 9. 

Elliptic cylinder - The flow around an elliptic cylinder at 45° incidence was 
calculated at Reynolds numbers of 200 and 3000. The computational mesh is 41x161 
for Reynolds number 200 and 81x241 for Reynolds number 3000. The far-field boundary 
is equivalent to 20 radii. The time evolution of the lift and drag coefficients is 
shown in figure 6 for Reynolds number 200. The time, t = 1.0, corresponds to the 
time required to travel a distance equal to the elliptic cylinder semi-major axis at 
a velocity equal to that of the free-stream velocity. The results compare favorably 
with the results of reference 10 up to t = 8. The exact trend in the data of 
reference 10 for t > 9.0 is unclear because those data were limited. The develop- 
ment in time of the streamlines around the elliptic cylinder are shown in figure 7 
for a Reynolds number of 200. The streamline contours are similar to those reported 
in reference 10. 

Variations in the lift and drag coefficients with time for the flow around an 
elliptic cylinder at Reynolds number 3000 are shown in figure 8. The cylinder is 
again at 45° incidence. The loading on the elliptic cylinder is shown to be more 
dynamic than that of the circular cylinder (see fig. 2). This is due to the large 
velocity gradient and resulting surface vorticity at the leading edge of the ellip- 
tic cylinder. In addition, the secondary vortex structure near the leading edge on 
the upper surface of the cylinder contributes significantly to the dynamic loading 
on the cylinder. The development in time of the streamlines around the cylinder are 
presented in figure 9. The development of the secondary vortices can be clearly 
seen. The separated wake behind the cylinder is shown to be very dynamic due to the 
strength and periodic shedding of vortices from the leading and trailing edges of 
the cylinder. 


Vortex Interaction Case 

The interaction of a vortex wake with the flow about an elliptic cylinder at 
45° incidence and a Reynolds number of 3000 was calculated. The computational mesh 
is 81x241 with the far field boundary at 20 radii. Two interactions were consid- 
ered, the difference between the interactions being the sign of the interacting 
vortex wake. Each interaction modeled the interacting vortex wake as an array of 
vortices inserted into the flow field at constant time intervals. This represents 
the periodic passage of a rotor blade and the associated tip vortex, thus generating 
the rotor's vortical wake structure. The orientation of the interacting vortex wake 
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with respect to the elliptic cylinder is depicted in figure 10. In each interaction 
the vortices are of constant strength and equal to r = -1.0 for the first inter- 
action and r = 1.0 for the second interaction. For both interactions r Q = 0.1. 
The flow was started impulsively at t = 0. Five equally spaced vortices are 
inserted into the flow field forward of the elliptic cylinder at t = 0.25. The 
flow around the elliptic cylinder is very unsteady; therefore, an optimum time to 
insert the vortices into the flow was not obvious. The time, t = 0.25, was chosen 
to minimize the effects of the impulsive start. A new vortex is inserted at a 
constant time interval, t = 1.0, during the remainder of the solution. A vortex is 
distributed to the mesh when the distance between the vortex position and the sur- 
face of the elliptic cylinder is less than 0.75. 

First Interaction, r = -1.0 - The lift and drag variations with time are shown 
in figures 11(a) and 11(b), respectively. Also shown on each figure is the respec- 
tive value of the lift and drag coefficient for the elliptic cylinder in the absence 
of the interacting vortex wake. Initially, for t < 0.25, the lift and drag coeffi- 
cients are identical to those calculated for the elliptic cylinder without the 
interacting vortex wake. However, after the interacting vortex wake is inserted 
into the flow, t = 0.25, the loading on the cylinder changes drastically. 

These changes seem to occur at t = 0.5 with a gradual decrease in the lift, 
figure 11(a), and drag, figure 11(b). At this time, the interacting vortex wake is 
below the cylinder and does not significantly affect the flow near the cylinder 
which is dominated by vortices generated by the cylinder during the impulsive 
start. As the interaction proceeds for t s 1.0, the lift (fig. 11(a)) and the drag 
(fig. 11(b)) are significantly reduced from those in the absence of the interacting 
vortex wake. The signature of individual vortices comprising the interacting vortex 
wake can be seen as spikes in the development of the lift coefficient (fig. 11(a)) 
at times t = 2.0, 3.0, 5.0, 7.0, and 8.5. Similar spikes can be seen in the devel- 
opment of the drag coefficient (fig. 11(b)) although less pronounced. 

The streamlines during the interaction are shown in figure 12. Initially, 
figure 12(a), the flow field is similar to that in the absence of the interacting 
vortex wake (fig. 9(a)) with the exception that the incoming streamlines have been 
rotated clockwise and the forward stagnation point has moved toward the leading 
edge. These flow changes are caused by the induced velocity of the interacting 
vortex wake and result in a reduction in the effective angle of attack of the ellip- 
tic cylinder. A reduction in the angle of attack decreases the strength of the 
separated vortices at the leading edge of the cylinder and changes the position of 
the secondary vortices. The secondary vortex structure is positioned farther down- 
stream of the leading edge (figs. 12(b-d)) when compared to the flow with no vortex 
interaction (figs. 9(b-d)). 

At later times (figs. 12(f-i)), the significant aspect of the flow is the 
generation of a strong vortex at the trailing edge of the cylinder. This vortex is 
caused by the increased flow around the trailing edge from lower to upper surface 
due to the interacting vortex wake and results in a very dynamic separated wake 
behind the cylinder. 
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Second interaction, r = 1.0 - The lift and drag variations with time shown in 
figures 13(a) and 13(b), respectively. Here again, the respective value of the lift 
and drag coefficient for the elliptic cylinder in the absence of the interacting 
vortex wake are shown for comparison. Initially, for t < 0.25 the loading on the 
cylinder is identical to that for the elliptic cylinder without the interacting 
vortex wake. Initial changes in the loading on the cylinder during the interaction 
from that with no vortex interaction occur at t = 0.5. 

These changes begin with a rapid increase in the lift (fig. 13(a)) and drag 
(fig. 13(b)). At this time, the interacting vortex wake is upstream of the cylin- 
der. However, the induced velocity caused by the interacting vortices generates an 
upwash at the leading edge of the cylinder. This upwash increases the relative 
angle of attack of the cylinder, and therefore the lift and drag on the cylinder. 

As the interacting vortex wake approaches the cylinder (t = 1.0) the lift 
(fig. 13(a)) and the drag (fig. 13(b)) are significantly increased from those in the 
absence of the interacting vortex wake. The signature of individual vortices com- 
prising the interacting vortex wake can be seen as spikes in the development of the 
lift coefficient (fig. 13(a)) at times t = 1.5, 2.5, 3.5, 4.8, 6.5, 7.6, and 9.6. 
Similar spikes can be seen in the development of the drag coefficient (fig. 13(b)). 

The streamlines produced during the interaction are shown in figure 14. Ini- 
tially (fig. 14(a)), the flow field is similar to that in the absence of the inter- 
acting vortex wake (fig. 9(a)) and of the previous interaction (fig. 12(a)); how- 
ever, the incoming streamlines have been rotated counterclockwise thus increasing 
the effective angle of attack of the cylinder. The increased angle of attack moves 
the forward stagnation point toward the midchord position and increases the flow 
around the leading edge of the cylinder. The increased flow at the leading edge of 
the cylinder increases the strength of the separated vortices generated at that 
position. In addition, the secondary vortex structure is positioned nearer the 
leading edge of the cylinder (figs. I4(b-d)) when compared to the flow with no 
vortex interaction (figs. 9(b-d)) and to the flow of the previous interaction 
(figs. 12(b-d) ) . 

As the interaction proceeds (figs. 1 4( f — i ) ) , the generation of the separated 
vortices at the leading edge of the cylinder dominates the flow. 


CONCLUSIONS 


A method has been developed to model the two-dimensional interaction between an 
interacting vortex wake represented by a finite-core model and the viscous flow 
around arbitrary bodies. The method solves for the flow field velocities on a body- 
fitted computational mesh using finite-difference techniques. The viscous flow 
field of the two-dimensional body is calculated on an Eulerian grid via the 
velocity/vorticity formulation of the Navier-Stokes equations. 

A simulation of a rotor wake interaction with the flow about a 4:1 elliptic 
cylinder at 45° incidence was shown for Reynolds number 3000. The simulation 
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demonstrates the significant variations in the lift and drag on the elliptic cylin- 
der in the presence of the interacting rotor wake. The interacting vortex wake 
significantly changes the effective incidence of the elliptic cylinder, thereby 
altering the loading on the elliptic cylinder and the generation of vorticity at the 
surface of the cylinder. 

When the strength of the interacting vortex wake was less than zero (producing 
a clockwise rotating vortex), the effective incidence of the cylinder (and therefore 
the loading) was significantly reduced from that for the flow around the cylinder in 
the absence of the interacting vortex wake. An interacting vortex wake of strength 
greater than zero had the opposite effect on the flow around the cylinder. The 
incidence of the cylinder and the loading on the cylinder were greatly increased. 

During both interactions, several secondary vortices were generated just aft of 
the leading edge on the upper surface of the cylinder. Oscillations in the loading 
on the elliptic cylinder during both interactions were shown. These oscillations 
were not present in the calculations of the flow in the absence of the interacting 
vortex wake and are felt to be the result of the increased shedding of secondary 
vortices. 
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Figure 1.- Computational geometry. 
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Figure 4.- Velocity profile downstream of circular cylinder, comparison with 
experiment, u velocity component, Re = 3000. 
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Figure 7.- Streamline contours, elliptic cylinder, Re = 200, 0 = 45 
(a) t = 1.00, (b) t = 2.00, (c) t = 3-00, (d) t = 4.00, (e) t = 5 
(f) t = 6.00, (g) t = 7.00, (h) t = 8.00, (i) t = 9.00, and (j) t 
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Figure 10.- Schematic of vortex interaction. 
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Figure 11.- Vortex interaction, vortex strength r = -1.0, r Q = 0.1, elliptic 
cylinder, Re = 3000, e s 45°. (a) Lift coefficient; and (b) drag coefficient 





Figure 12.- Streamline contours for a vortex interaction, vortex strength 
r s 1.0, r Q * 0.1, elliptic cylinder, Re = 3000, e = 45°. (a) t = 1.00 

(b) t = 2.00; (c) t = 3.00; (d) t = 4.00; (e) t = 5.00, (f) t = 6.00, 
(g) t s 7.00, (h) t = 8.00, (i) t = 9.00, and (j) t = 10.00. 
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VORTEX INTERACTION 



Figure 13.- Vortex interaction, vortex strength r = 1 .0, r Q = 0. 1 , elliptic 
cylinder, Re = 3000, 0 = 45°. (a) Lift coefficient; and (b) drag coefficient. 
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